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The average-atom model is applied to study Thomson scattering of x-rays from warm-dense mat- 
ter with emphasis on scattering by bound electrons. Parameters needed to evaluate the dynamic 
structure function (chemical potential, average ionic charge, free electron density, bound and con- 
tinuum wave functions and occupation numbers) are obtained from the average-atom model. The 
resulting analysis provides a relatively simple diagnostic for use in connection with x-ray scatter- 
ing measurements. Applications are given to dense hydrogen, beryllium, aluminum and titanium 
plasmas. In the case of titanium, bound states are predicted to modify the spectrum significantly. 
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I. INTRODUCTION 

Measurements of Thomson scattering of x-rays pro- 
vide information on temperatures, densities and ioniza- 
tion balance in warm dense matter. Various techniques 
for inferring plasma properties from x-ray scattering mea- 
surements have been developed over the past decade [H- 
|22| ; these techniques together with the underlying theory 
were reviewed by Glenzer and Redmer [23[ 

The present analysis of Thomson scattering from dense 
plasmas is based on a theoretical model proposed by Gre- 
gori et al. Q , one important difference being that the pa- 
rameters used here to evaluate the Thomson-scattering 
dynamic structure function are taken from the average- 
atom model. The particular average-atom model used 



here is described in Ref. [24( . The present work is closely 
related to that of Sahoo ct al. [25(, where a somewhat 
different version of the average-atom model was used. 
Predictions from t he p resent model differ substantially 
from those in Ref. [25[. The origin and consequences of 
these differences will be discussed later. 

The Thomson scattering cross section for an incident 
photon with energy, momentum (ujq, fco) and polariza- 
tion eo scattering to a state with energy, momentum 
(wi, fei) and polarization e\ is 
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The dynamic structure function S(k,ui) appearing in 
Eq. (fT|) depends on two variables: k = \ko — and 
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ui = luq — uj\. As shown in the seminal work of Chi- 
hara [2(| [27j, S(k,ui) can be decomposed into three 
terms: the first Sa(k,uj) is the contribution from elas- 
tic scattering by electrons that follow the ion motion, 
the second S ee {k,uj) is the contribution from scattering 
by free electrons and the third SB(k,uj) is the contri- 
bution from bound-free transitions (inelastic scattering 
by bound electrons) modulated by the ionic motion. In 
the present work, the modulation factor is ignored when 
evaluating the bound-free scattering structure function. 
For the bound-free contribution, calculations carried out 
using plane-wave final states are compared with calcu- 
lations carried out using average-atom scattering wave 
functions. Substantial differences are found between 
these cases. 

The average-atom model is discussed briefly in Sec. UII 
followed by a discussion of the three contributions to the 
structure functions in Scc. lIIII In Scc. lIVI applications are 
given to hydrogen, beryllium, aluminum and titanium 
plasmas. 



II. AVERAGE-ATOM MODEL 

Average-atom models arc versions of the temperature- 
dependent Thomas-Fermi model of a plasma developed 
63 years ago by Feynman et al. [28| which include de- 
tailed descriptions of bound and continuum states of an 
atom imbedded in a plasma. In this model, the plasma is 
divided into neutral Wigner-Seitz (WS) cells (volume per 
atom V ws = A/pNa, where A is the atomic weight, p is 
the mass density and Na is Avogadro's number). Inside 
each WS cell is a nucleus of charge Z and Z electrons. 
Some of these electrons are in bound states and some 
in continuum states. The continuum density is finite 
at the cell boundary and merges into the uniform free- 
electron density Zf/V ws outside the cell. Each neutral 
cell can, therefore, be regarded as an ion imbedded in a 
uniform sea of free electrons of density n e = Zf/V ws . To 
maintain overall neutrality, it is necessary to introduce a 
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uniform (but inert) positive background density Zf/V ws . 
The model, therefore, describes an isolated (neutral) ion 
floating in a (neutral) "jcllium" sea. 

The quantum-mechanical model here, which is dis- 
cussed in Ref . j24| , is a nonrelativistic version of the rel- 
ativistic Inferno model of Liberman [29fl and the more 
recent Purgatorio model of Wilson et al. [3 Of ; it is similar 
to the nonrelativistic average-atom model described by 
Blenski and Ishikawa [3l[. Specifically, each electron in 
the ion is assumed to satisfy the central-field Schrodinger 
equation 



Z 
r 



V 



-0a( r ) = e a 1pa{r), 



(3) 



where a = (n, I) for bound states or (e, I) for continuum 
states. Atomic units (a.u.) where e = % = m = 47reo = 

1 are used here. In particular, 1 a.u. in energy equals 

2 Rydbergs (27.211 eV), and 1 a.u. in length equals 1 
Bohr radius a (0.529 A). The wave function tp a (r) is 
decomposed in a spherical basis as 



Mr) = -Pa(r)Y lama (f) Xaa , 



(4) 



where Yi m (f) is a spherical harmonic and Xn is a two- 
component electron spinor. The bound and continuum 
radial functions P a (r) are normalized as 



drPniir) P n >i(r) = 5 r , 



drP d (r) P en {r) = S(e 



(5) 
(6) 



respectively. The central potential V(r) in Eq. (HI is 
taken to be the self-consistent Kohn-Sham potential [32[ 



V(r) = 4tt / — r' 2 n(r') dr' - x a 



81 . 
-n(r) 



(7) 



where the first term in the right-hand side of Eq. (|7|) is 
the direct screening potential with r> = max(r, r') and 
the second term is the average exchange potential with 
x a = 2/3. While short-range electron-electron interac- 
tions inside the Wigncr-Scitz cells are reasonably well 
accounted for by this simple model, it should be noted 
that eigenvalues in the Kohn-Sham potential are poor 
approximations to ionization energies, leading to inaccu- 
rate thresholds and peaks of bound-free contributions to 
S(k,u>), which can be off by 20 - 30% when compared 
with experiment. 

The electron density n(r) in Eq. ([7]) has contribu- 
tions from bound-states n b (r) and from continuum states 
n c (r), 

n(r) = n b (r) + n c (r). (8) 
The bound-state contribution to the density n b {r) is 

2(2/ + 1) 

1 4- ex 
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+ exp[(e„; - fi)/k B T] 



Pm{r) 2 , (9) 
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FIG. 1: (Color online) Upper panel: The radial density 
4nr 2 n(r) for Al at metallic density and k B T — 5 eV (solid 
curve) integrates to Z = 13 for r < Rws- The contin- 
uum contribution 47rr 2 n c (r) (dashed curve) integrates to 3 
for r < i?ws- The bound Is, 2s and 2p shells are completely 
occupied at this temperature. Lower panel: The dashed curve 
illustrates the Friedel oscillations of the continuum density 
and shows how N c (r) — n c (r)Vws converges to Zf = n e Vws 
(solid line) for r > i?ws- The chemical potential predicted by 
the model is /j, = 0.2406 a.u. and the number of free electrons 
per ion is Zf = 2.146. 



where e„; is the bound-state energy, \x is the chemical po- 
tential, and the sum over (n, I) ranges over all bound sub- 
shells. The continuum contribution to the density n c (r) 
is given by a similar expression with the bound state 
radial functions P n i{r) replaced by continuum functions 
Pei{r) and the sum over n replaced by an integral over 
e. Finally, the chemical potential fj, is chosen to ensure 
charge neutrality inside the WS cell: 



Z= I n(r)d"r 

'r<i? w s 



3„ — 



-Rws 



4wr 2 n(r) dr 



(10) 



Equations |3l fT0|) above are solved self-consistcntly to give 
the chemical potential /i, the potential V(r) and the elec- 
tron density n(r). Numerical details can be found in 
Ref. [13. 

The boundary conditions used in solving Eq. ([3]) de- 
serve some mention. Bound state wave functions and 
their derivatives are matched at the boundary r = i? ws 
to solutions outside the WS sphere (where V = 0) that 
vanish exponentially as r — > oo. Similarly, continuum 
functions and their derivatives are matched to phase- 
shifted free-particle wave functions at r = i? W s ■ It should 
be noticed that the continuum density n c (r) inside the 
WS sphere, which oscillates as predicted by Friedel (33| . 
is distinctly different from the uniform free electron den- 
sity n e . In the present model, n c (r) smoothly approaches 
n e outside the sphere. These points are illustrated in 
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Fig. Q] where the bound-state and continuum densities 
are plotted for Al at metallic density and temperature 
k B T = 5 cV. 

The boundary conditions used here differ from those 
used by Sahoo et al. in Ref. (25|, where the first deriva- 
tive of the wave function is required to vanish at i? W s- 
The differences in boundary conditions lead to major 
differences in the average-atom structure. For example, 
the model used in [25[ predicts that the M shell of Al 
at metallic density is partially occupied at temperatures 
k B T < 10 eV, whereas the present model predicts that 
the M shell is empty in this temperature range. Conse- 
quences of such differences are discussed later in Sec. IIVI 

III. DYNAMIC STRUCTURE FUNCTION 

In the paragraphs below, the evaluation of S(k,w) in 
the average-atom approximation is discussed. As men- 
tioned earlier, the theoretical model developed by Gre- 
gori et al. Q is used to evaluate the ion-ion contribu- 
tion Sa(k,ui) to the dynamic structure function. Addi- 
tionally, the procedure proposed in Ref. @ is used to 
account for differences between electron and ion tem- 
peratures. The electron-electron contribution S ee {k,ui) 
is expressed in terms of the dielectric function e(fc,w) 
of the free electrons which in turn is evaluated using 
the random-phase approximation (RPA) . Finally, bound- 
state contributions to the dynamic structure function are 
evaluated using average-atom bound state wave func- 
tions. The final-state wave function is described in two 
different ways: (1) using a plane-wave final-state wave 
function as in Ref. [25[, and (2) using an average-atom 
final-state wave function that approaches a plane wave 
asymptotically. There are dramatic differences between 
these choices. The more realistic average-atom choice 
automatically includes ionic Coulomb-field effects. 

A. Ion-Ion Structure Function 

The contribution to the dynamic structure function 
from elastic scattering by electrons following the ion mo- 
tion Sa{k,u>) is expressed in terms of the corresponding 
static ion-ion structure function Su{k) as: 

Sa{k,u) = |/(fc) + <z(fc)| 2 S ti (k)6(u). (11) 

In the above, f(k) is the Fourier transform of the bound- 
state density and q(k) is the Fourier transform of elec- 
trons that screen the ionic charge. In the average-atom 
approximation, the screening electrons are the continuum 
electrons inside the Wigner-Seitz sphere and 

rRws 

f(k) + q(k) = 4ir r 2 [n b {r)+n c (r)]j (kr)dr, (12) 

Jo 

where ji(z) are spherical Bcsscl functions of order I. Note 
that /(0) + q(0) = Z in the average-atom model. Fur- 
thermore, the delta function S(lu) in Eq. (jll[) is replaced 
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FIG. 2: (Color online) Upper panel: Sa(k) is shown for Be 
metal at electron temperature T e =20 eV and ion-electron 
temperature ratios Ti/T e =(l, 0.5, 0.1) illustrated in solid, 
short dashed and long dashed curves, respectively. The value 
k — 0.543 corresponds to an incident photon uia = 2960 eV 
scattered at angle 40°. Lower panel: Sa(k,u) for Be metal 
at T e = 20 eV and Tj = 2 eV, where the function S(lj) is 
replaced by a Gaussian of width 10 eV and k — 0.543. 



by an "instrumental" Gaussian, with full-width at half 
maximum = 10 eV in this work. This value is chosen be- 
cause typical experiments in Be [l3| have a spectrometer 
with a 10 eV instrument width and use a CI Ly-cv source 
at 2.96 keV. 

Approximate schemes to evaluate the static structure 
functions Sa(k) are discussed, for example, in Ref. [34l ]. 
Here, we follow Ref. [|| and make use of formulas given 
by Arkhipov and Davletov [3{| that account for both 
quantum-mechanical and screening effects. The func- 
tion Su(k) in Ref. (3f| is expressed in terms of the 
Fourier transform of the ion-ion interaction potential 
<&ii(r) through the relation: 

77 - 

S«(*0 = 1 " -Hr*«(*). ( 13 ) 
k B T 

where n.j is the ion density. 

Different Electron and Ion Temperatures In the aver- 
age atom model, T is the electron temperature T e which, 
in equilibrium, is equal to the ion temperature X^. To al- 
low for different electron and ion temperatures, the equa- 
tions for Su(k) given by Arkhipov and Davletov (3oT | are 
modified following the prescription laid out by Gregori 
et al. [H . The electron temperature T e is replaced by an 
effective temperature T' e that accounts for degeneracy ef- 
fects at temperatures lower than the Fermi temperature 
Tp. Similarly, the ion temperature Tj is replaced by an 
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FIG. 3: (Color online) Upper panel: Real (solid) and 
imaginary (dashed) parts of e{k,ui) are plotted along with 
Im[l/e(fc,w)] (dot-dashed) for Be metal at k B T = 20 eV. 
Lower panel: The resulting structure function S ee (k,io) 
(solid) is shown together with Su(k,u)) (dashed). These plots 
correspond to Thomson scattering of a 2960eV photon at 20° . 



effective temperature T[ that accounts for ion degener- 
acy effects at temperatures lower than the ion screened 
Debye temperature Trj. Explicit formulas for Su(k) are 
found in Ref. [|| . The dramatic effect of different electron 
and ion temperatures on the static structure functions 
Su(k) for Be at metallic density and T e = 20 eV are il- 
lustrated in the top panel of Fig. [5J This figure is similar 
to the upper-left panel of Fig. 1 in Ref. which was 
obtained under similar condition. In the bottom panel 
of Fig. [2] contributions to Su{k, uj) for Be at T e = 20 eV 
and T; = 2 eV are shown. 



B. Electron-Electron Structure Function 

The electron-electron structure function S ee {k,ui) is 
expressed in terms of the plasma dielectric function 
e(k,uj) through Eq. (15) in Ref. @: 



S ee (k,uj) 



1 



1 — exp(— ui/k B T) 4irn e 



Im 



1 



e{k,u) 



(14) 

In the average atom model, the free electrons are uni- 
formly distributed outside the WS sphere. The density of 
these electrons is n e = Zf /V ws . In the present work, the 
dielectric function is evaluated using the random-phase 



approximation. The real and imaginary parts of the RPA 
dielectric function e(fc,w), given in Eq. (16) of Ref. H, 
can be written as 



Rc[e(fc,w)] 



In 
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(15) 



and 



2 f b 

Im[e(k, uj)] = -3 / F(p) p dp 

-exp[(fi~a 2 /2)/k B T] 
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with a = \2ui 
equations, 



k 2 \/2k and b = (2u> + k 2 )/2k. In these 



Hp) 



l + exp[(p 2 /2-/i)/fc B T] 



(17) 



is the free-electron Fermi distribution function. It should 
be noted that the real part of e(k,uj) is an even function 
of ui and the imaginary part is an odd function of uj. 

The real and imaginary parts of e(k,uj) along with 
— Im[l/e(fc,o;)] are illustrated in the top panel of Fig. [3] 
for scattering of a 2960 eV photon at 20° from Be metal 
at 20 eV. The sharp peak in Im[l/e] that occurs near 
the point where Re[e] vanishes is a collective plasma 
resonance (plasmon). The contribution to S ee (k,ui) is 
shown in the bottom panel. The ratio of the down-shifted 
(u>i < 2960 eV) to up-shifted (uii > 2960 eV) resonance 
peaks exp(Aw/fc B T), where Auj is the energy of the plas- 
mon peak relative to the central energy, is used to deter- 
mine the electron temperature. 



C. Scattering from Bound States 

The structure function associated with Thomson scat- 
tering from a bound state ipniir) with quantum numbers 
(n, /) to a continuum state ipp{ r ) with momentum p is 



S„i(k,u>) 



E 



pdD, p 



(2tt) 2 



d\^ p {r)e lk - r ^ nlm (r) 

. (18) 

As mentioned previously, two possibilities are considered 
for the final state in bound-free scattering: (1) a free- 
particle plane wave, and (2) an average-atom contin- 
uum wave that approaches a plane wave asymptotically. 
Case (2) is clearly the more physical alternative since 
continuum waves in the average-atom potential differ 
markedly from free-particle wave functions. This point 
is illustrated in Fig. [U where the average atom radial- 
functions P e i (r) jpr are compared with their free-particle 
counterparts ji(pr). The average-atom wave functions 
are seen to differ markedly from the free-particle (spher- 
ical Bessel) functions for low values of I, but approach 
free-particle functions as I increases. 
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FIG. 4: (Color online) Comparison of average-atom contin- 
uum functions Pi(pr)/pr (solid lines) for Al at metallic den- 
sity and k B T — 5 eV with spherical Bessel functions ji(pr) 
(dashed lines). 
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FIG. 5: (Color online) The beryllium K shell structure func- 
tion Sis (k, uS) is shown for incident photon energy 2960 eV and 
scattering angles 30° and 150°. The black curves show plane- 
wave results, the blue lines show the result obtained using an 
average atom final-state wave function and the red lines show 
exact nonrelativistic Coulomb results. The dramatic suppres- 
sion of average atom and Coulomb structure functions at for- 
ward angles (the corresponding curves are multiplied by 10) 
is evident in the upper panel. 



a. Plane-wave final states Assuming that the final 
state wave function is a free-particle plane wave e lp r , the 
bound-free structure function in Eq. (|18p can be rewrit- 
ten as 



S n i(k,uj) = 



p cM p 



d''i 



(19) 

where fe — fco — fei, w = loo — and q = k p. Note 
that q is the momentum transferred to the ion. This 
expression may be simplified to 



Sm{k,u) = / qdq\K nl (q)\ 2 , (20) 

-l\p-k\ 



p+k 



where o n i is the occupation number of the final state and 

/>oo 

K n i(q) = / drrji(qr)P n i(r). (21) 



Eq. (|20[) depends implicitly on uj through the relation 



p = \/2(uj + e n i) . 

b. Average-atom final states In the average-atom 
approach, the final state wave function consists of a plane 
wave plus an incoming spherical wave. (n.b. An outgo- 
ing spherical wave is associated with an incident electron. 
Time-reversal invariance, therefore, requires that a con- 
verging spherical wave be associated with an emerging 
electron.) The bound- free structure function in Eq. (|18[) 



may be reexpressed as 

S n i = -0^411, \I h ii 2 (p,k)\ 2 , (22) 



hi. 



where o n i is the occupation number of the final state with 

n-RwS 



1 rftWS 

Ihih (P, k)=- e*^' / air P dl (r) Jh (kr)P nl (r). 
P Jo 



and 



(23) 



A llll2 = {2h + l){2h + l)^ I (24) 

In the above, Si ± (p) is the phase-shift of the final state 
partial wave (r). Moreover, e = ui + e n i, p = v/2e and 
k = \k Q - fei|. 

In Fig. [5l several calculations of the structure function 
S\ s (k,uj) are compared for a photon of incident energy 
2960 eV scattered at 30° and 150° from the K shell of 
beryllium metal at T = 20 eV. The results of calculations 
carried out using average-atom final states are smaller 
than those using plane-wave final states by a factor of 
about 40 at forward angles and 2.5 at backward angles. 
This suppression is a characteristic Coulomb field effect. 
Indeed, exact nonrelativistic Coulomb-field calculations 
of Thomson scattering |36| . with nuclear charge adjusted 
to align the Coulomb and average atom thresholds, show 
a similar suppression. 
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FIG. 6: (Color online) The evolution of plasmon resonances 
for scattering of a 5 keV photon from a fully ionized dense H 
plasma (T = 50 eV, n e = 10 24 cm" 3 ) is illustrated in the top 
panel where we plot S(k,uj) for scattering angles of 20°, 30° 
and 40°. The corresponding coherence parameters a — 1/X s k 
are plotted in the bottom panel. 



IV. APPLICATIONS 

In the subsections below, S(k,Lj) is evaluated in the 
average-atom approximation for cases of possible ex- 
perimental interest: hydrogen at n e = 10 24 cm -3 and 
T = 50 eV, beryllium at n e = 1.8 x 10 23 cm -3 and 
T = 18 eV, aluminum at metallic density and T — 5 eV 
and titanium at metallic density and T = 10 eV. 



A. Hydrogen at T = 50 eV and n e 



10 2 



In the average-atom model, a density p =1.931 g/cc is 
required at T = 50 eV to achieve free-electron density 
n e = 10 24 cm -3 . The chemical potential in this case is 
p = —1.091 (a.u.). Under these conditions of temper- 
ature and density, hydrogen is completely ionized. The 
continuum density n c (r) inside the WS sphere merges 
into the free-electron density n e outside the sphere. The 
total number of electrons inside the WS sphere N c — 
An J flwb r 2 n c (r)dr = 1, however, Zf = 0.8667. 

Since there are no bound electrons in this case, only 
Su and S ee contribute to the cross section. Dynamic 
structure functions for scattering of a 5 keV photon at 
angles ranging from 20°, 30° and 40° are shown in the 
top panel of Fig. [5] Resonance peaks are seen to broaden 
and move to higher frequencies as the scattering angle 
increases. The coherence parameter a = l/(A s fc), defined 



in Eqs. (5-7) of Rcf. [23|], is plotted in the bottom panel 
of Fig. [5] The parameter A s is the shielding length, given 
by 



/ k B TF 1/2 (p/k B T) 
3 Y ^n e F_ 1/2 (p/k B T)' 

where Fj(x) is a complete Fermi-Dirac integral, 



(25) 



(26) 



For this particular case, A s = 1.071 a.u.. The value of 
A s differs only slightly from the WS radius R ws = 1.118 
a.u.. The resonant features in Fig.|5]are distinct for a > 1 
but disappear for a < 1, in harmony with the fact that 
plasmon resonances are collective phenomena. It should 
be noted that the (unperturbed) plasma frequency for 



hydrogen at n e = 10 cm 



j p i = 37.1 eV. 



B. Beryllium at T = 18 eV and n e = 1.8 xlO 23 cm -3 

In the bottom panel of Fig. [71 the structure function 
for scattering of a 2963 eV photon at 40° from beryl- 
lium (density = 1.636 g/cc) at T e = 18 eV is plotted. 
The L shell electrons are completely stripped under these 
conditions but the K shell remains 97% occupied. The 
chemical potential is p = —0.5311 a.u. and the number 
of free electrons per ion Zf = 1.647. The ion tempera- 
ture, which governs the amplitude of the elastic peak, is 
chosen to be T, = 2.1 eV in this example. For the case 
at hand, the coherence parameter is a = 1.21, so one 
expects and observes plasmon peaks in the scattering in- 
tensity profile. The average-atom removal energy for a 
K shell electron is 86.8 eV. One therefore expects to find 
a contribution to S(k,u>) from K shell electrons for ener- 
gies uji < 2876 eV. The K shell contribution multiplied 
by 50 is shown in the bottom panel. 

To validate the present average-atom model against 
experimental data, a Be experiment done at the Omega 
laser facility that used a CI Ly-a source to scatter from 
nearly solid Be at an angle of 40° is used. An electron 
temperature of 18 eV, ion temperature of 2.1 eV, and 
density of 1.647 g/cc used in the average-atom model 
gives an electron density of 1.8 x 10 23 /cc, in agreement 
with the analysis in Ref. [13 • The top panel of Fig. [7] 
shows the experimental source function from the CI Ly- 
a line as a blue dashed line. Because of satellite structure 
in the source we approximate the source by three lines: 
a CI Ly-a line at 2963 eV with amplitude 1 and two 
satellites at 2934 and 2946 eV with relative amplitudes 
of 0.075 and 0.037 respectively. Doing the Thomson scat- 
tering calculation using the three weighted lines, we cal- 
culate the scattering amplitude for Thomson scattering 
(red dashed line) and compare against the experimen- 
tal data (black solid line) here. We observe excellent 
agreement within the experimental noise. Contributions 
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FIG. 7: (Color online) Bottom panel: Structure function 
S(k, oj) for scattering of a 2963 eV photon at 40° from beryl- 
lium at metallic density and T e = 18 eV. Top panel: Intensity 
measurements for scattering of a CI Ly-a x-ray from beryl- 
lium at 40° (37|; measurement (solid line), source function 
(dot-dashed line), and average atom fit (dashed line). 



from the bound Is electrons, which have a threshold at 
2876 eV, are beyond the range of the data shown in the 
top panel. 



C. Aluminum at T = 5 eV and metallic density 

Aluminum at metallic density (p = 2.70 gm/cc) and 
T = 5 eV has a Nc-like ion configuration with two 2s 
electrons bound by 92.2 eV and six 2p electrons bound by 
54.9 eV. There are three continuum electrons inside the 
WS sphere i? W s = 2.99 (a.u.). The continuum density 
inside the sphere n c (r) converges to the uniform free- 
electron density n e = Zf/V WK , where Zf = 2.146. In 
Fig. HI the structure function S(k,u) is plotted for the 
case of an incident 2.96 keV photon scattered at 30°. The 
coherence parameter a = 1.95 in this case, explaining the 
prominent plasmon resonance seen on the low-frequency 
side of the elastic scattering peak. Also shown in the 
figure are contributions from the bound L shell electrons 
scaled up by a factor of 100. It should be noted that, 
by contrast with the average-atom calculations presented 
in Ref. [25[, the M shell of Al is completely empty at 
temperatures below 10 eV in the present model and the 
prominent M shell features predicted in Ref. [25| do not 
arise in the present analysis. 



FIG. 8: (Color online) Structure function S(k, u) for scat- 
tering of a 2960 eV photon at 30° from metallic density alu- 
minum (Z = 13) at T e = 5 eV and Ti/T e = 0.1. Contribu- 
tions from the L shell of the Ne-like core multiplied by 100 
are indicated on the plot. 



D. Titanium at T — 10 eV and metallic density 

Titanium (Z = 22) at metallic density (p = 4.51 g/cc) 
and T = 10 eV is a case where sharp peaks from excita- 
tions of bound M shell electrons show up in the energy 
spectrum along with the plasmon peaks. The average- 
atom model predicts that metallic density Ti is in an 
Ar-likc configuration at T = 10 eV with completely filled 
K and L shells together with 1.97 3s electrons bound 
by 44.40 eV and 5.36 3p electrons bound by 22.88 eV. 
There are 4.67 continuum electrons inside the WS sphere 
i? ws = 3.05 (a.u.). The continuum density inside the 
sphere converges to the uniform free-electron density 
n e = Zf/V WK outside the sphere, with Zj = 2.305. The 
chemical potential is p = —0.0511 au. In Fig. [SJ the dy- 
namic structure function S(k, ui) is shown for the case of 
an incident 2.96 keV photon scattered at 30° and 150°. 
Plasmon peaks, which are prominent for scattering at 
30°, disappear for scattering at 150° while the 3s and 
3p bound-state peaks grow. The M shell contributions 
to the structure function are comparable to the plasmon 
contribution for the 30° case and are the dominant fea- 
tures on the low-frequency side of the elastic peak at 
150°. 



V. SUMMARY 

A scheme for analysis of Thomson scattering from 
plasmas based on the average-atom model, a quantum- 
mechanical version of the "Generalized Thomas-Fermi 
Theory" of Feynman, Metropolis and Teller [28j is pre- 
sented. Given the plasma composition (Z,A), density p 
and temperature T, the model gives, in addition to the 
equation of state of the plasma, all parameters needed for 
a complete description of the Thomson scattering pro- 
cess. In particular, the average-atom code predicts wave 
functions for bound and continuum electrons, densities 
of bound, screening and free electrons, and the chemical 
potential. 
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FIG. 9: (Color online) Dynamic structure function S(k, u) 
for scattering of a 2960 eV photon at 30° and 150° from 
metallic density Ti at T e = 10 eV. Contributions to S(k,ui) 
(solid curve) from Se e (k, to) (dashed curve) and 5b (fc, u) (dot- 
dashed curve) are shown. 



is treated following the model proposed by Gregori et al. 
Q which makes use of formulas for the static ion-ion 
structure function Sa(k) given by Arkhipov and Davle- 
tov [35|. Modifications suggested by Gregori et al. j|| 
to account for different electron and ion temperatures 
are also included. Treatment of the ion-ion structure 
function appears to be the weakest aspect of the present 
analysis. The dynamic structure function for scatter- 
ing from free electrons depends sensitively on the free- 
electron dielectric function e(fc,w). Again, we follow the 
model proposed in Ref. Q and evaluate the dielectric 
function in the random-phase approximation. The RPA 
dielectric function includes features such as plasmon res- 
onant peaks that show up in experimental intensity pro- 
files and can be used in connection with the principle 
of detailed balance to determine electron temperatures. 
Bound-state features are included in the present scheme, 
inasmuch as the average-atom model provides bound- 
state and continuum wave functions. Coulomb-field ef- 
fects are automatically included in calculations carried 
out using average-atom continuum states rather than 
plane waves to describe the final state electron. In con- 
clusion, the average-atom model provides a simple and 
consistent point of departure for theoretical analysis of 
Thomson scattering from plasmas. 



Predictions of the present average-atom model disagree 
with those in Ref. [25[ where a similar model with dif- 
ferent boundary conditions is used. In particular, in 
Ref. [25|, 3d electrons were bound in metallic density Al 
for temperatures between 2 and 10 eV, leading to sub- 
stantial bound-state contributions to the dynamic struc- 
ture function. In the present model the 3d subshell of 
metallic density Al is vacant in the temperature range 
T < 10 eV and the corresponding bound-state features 
are absent. 

Elastic scattering from bound and screening electrons 
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